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t I We produce a detailed map of the perturbed surface density of neutral hy- 

^ ■ drogen in the outer Milky Way disk demonstrating that the Galaxy is a non- 

^ '. axisymmetric multi-armed spiral. Spiral structure in the southern half of the 

^ I Galaxy can be traced out to at least 25 kpc, implying a minimum radius for 

^ ; the gas disk. Overdensities in the surface density are coincident with regions 

1 ■ of reduced gas thickness. The ratio of the surface density to the local median 

2 ■ surface density is relatively constant along an arm. Logarithmic spirals can be 

^ ' fit to the arms with pitch angles of 20°-25°. 

. 

>• . Mapping the Milky Way's spiral structure is traditionally difficult because the Sun is imbed- 

^ i ded in the Galactic disk; absorption by dust renders optical methods ineffective at distances 

i-j I larger than a few kpc. Radio lines like the 21 cm hyperfine transition of atomic hydrogen (HI) 

are not affected by this absorption, and are therefore well-suited to looking through the disk. 
The density of HI is roughly proportional to the intensity of the emission, barring optical depth 
effects. Maps are constructed by using the Doppler shift of the emission line in combination 
with the rotation structure to determine where in the Galaxy the emission originates. The first 
maps of the HI density in the midplane of the Milky Way { !) made using the 21 cm transi- 
tion offered the promise that HI mapping of the Galactic plane would reveal its spiral structure 
by highlighting regions with HI overdensities. The task is easiest outside of the Sun's orbit 
around the Galactic center because velocities in this region map into distances uniquely; in the 
inner Galaxy, however, each velocity corresponds to two distinct distances from the Sun, so an 
unambiguous velocity-to-position mapping is impossible. When observations of the southern 
plane of the Galaxy made in Australia were combined with northern sky data from the Nether- 
lands Q), a picture of spiral structure of the inner Milky Way was seen, but it did not clearly 
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resemble that of a typical spiral galaxy; the spiral structure in the outer Galaxy (outside the 
Solar orbit) was less apparent. More recent maps of HI in the outer Galaxy have improved on 
the early work {3). 

In this paper we present a map of perturbations on the HI surface density in the outer Milky 
Way using a modified unsharp masking technique. By subtracting a blurred copy from the 
original image, this technique emphasizes low contrast features in both bright and dim regions 
(0) such as spiral structure. In previous maps (e.g. (3) and Fig. 1) spiral arms were difficult to 
pick out on top of the surface density's global falloff with radius. 

We construct a grid of the HI density for the outer Galactic disk, p, as a function of Galacto- 
centric radius R, angle (p, and height off the plane z from the Leiden/ Argentine/Bonn Galactic 
HI survey (|5|0|3I11). The density is recovered from the measured intensity by assuming a con- 
stant spin temperature of 155 K and using the Doppler shift velocity to determine distance along 
the line of sight. We apply a median filter to remove small features not associated with the disk 
(such as other galaxies), and a thickness filter to remove larger-scale features like clouds and 
spurs coming off the disk. We adopt the International Astronomical Union standard value of 
Rq = 8.5 kpc for the Sun-Galactic center distance; recent measurements suggesting that this 
value may be too high would result in a proportional rescaling of the disk L9..10y We construct 
the surface density map, 0), by summing the grid p(R, 0, z) over z (Fig. 1). For our rota- 
tion structure, we use circular orbits with a frequency m = 1 epicyclic streamline correction and 
a circular velocity of 220 km/s constant with radius (JJJ. The epicyclic correction is necessary 
because previous studies of HI spiral arms that use purely circular orbits result in grossly dis- 
continuous surface density contours near the Sun-Galactic center line (e.g. ©); the magnitude 
of the correction is fit to minimize this discontinuity. The maximum line-of-sight velocities of 
the corrections are fixed to lie along the Sun-Galactic center line. We exclude from the map 
regions within ±15° of the Sun-Galactic center line because of the difficulty of determining 
reliable distances in these directions. Our radial grid runs over the Galactocentric range 8.6 kpc 
to 40 kpc. 

For each point on our grid, we calculate r]{R, 0), the median surface density of the points 
within 25° in and 2 kpc in R; this local median is analogous to the blurred image used in 
unsharp masking. It is well suited for looking for perturbations in surface density because it 
takes into account the falloff in surface density with radius, and adjusts if the Galaxy has a 
lopsided or lumpy distribution. For points within 2 kpc of the inner and outer radial borders of 
our grid, we narrow the radial range so that the median is evenly balanced with surface densities 
inside and outside the point. We then compare the surface density at each point to the median 
value: 

n(i?,0) = S(i?,0)/r/(i?,0). (1) 

This is a dimensionless quantity, and therefore is a direct measure of the strength of the sur- 
face density perturbations (Fig. 1). Dividing by the local median rather than subtracting, as is 
normally done in unsharp masking, compensates for the change in the surface density by more 
than an order of magnitude over the radial range of the map. 

Values of surface density ratios in Fig. 1 range from a minimum of 11 = 0.13 to a maximum 
of 10.36. The vast majority of points have values in the range 0.6 to 1.8, implying that the 
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Figure 1: A: Contour plot of the surface density S(_R, 0). The location of the Sun is marked 
at (0, 8.5 kpc) by the solar symbol 0. The regions excluded due to unreliable distances are the 
large blank wedges near the Sun-Galactic center line. B: Contour plot of Il{R, 0), as defined 
in the text. Colored regions are overdense compared to the local median, whereas grayscale 
regions are underdense. The solid contour marks the line 11 = 1.1. The values of 11 for the 
different contour levels are given by the colorbar. 
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arm-to-interarm surface density ratio is about 3. The typical value of 11 does not vary strongly 
as a function of R along an arm, or even from arm to arm. The 11 = 1 line does not signify 
the border of a spiral arm, so Fig. 1 does not provide a meaningful description of the width 
of the arms. Furthermore, the intrinsic velocity dispersion serves to widen the apparent arms, 
especially at large radii. 

Several long spiral arms appear clearly on the map, but the overall structure does not have 
the reflection through the origin symmetry of a "grand-design" spiral. Easily identified spiral 
arms cover a larger area in the south than in the north; there are three arms in the southern half of 
the diagram. Calculating r] over a smaller range of R and reveals more structure in the arms, 
but does not change their positions. Each of these arms has already been detected in previous 
maps out to about 17 kpc for the two arms closer to the anticenter (1) and out to 24 kpc for the 
'outer' arm (12). The modified unsharp masking technique allows us to trace the arms further; 
these arms run coherently over a length of nearly 30 kpc in the outer Galaxy alone. Near the 
Sun, the overdense region in Fig. 1 at i? ~ 10 kpc and ~ 160° is associated with the Perseus 
arm; however, there is a well-known discrepancy between photometrically and kinematically 
determined distances in this region ai3\U4i . The well-defined underdense interarm regions in 
the south are also quite striking. There are clear underdense regions between all three of the 
arms in the south. In the north, a strongly underdense region is apparent inside the arm at 
R^13 kpc. 

When using the Doppler shift velocity of an emission line to find the distance to the emitting 
gas, small-scale velocity perturbations result in distortions of the distance measure because the 
assumed rotation structure is inaccurate. In fact, peaks in 21 cm emission intensity are more 
likely to result from a combination of density and velocity perturbations (17511731 than solely 
from HI overdensities. Physically, this occurs because a peak can be explained by a local density 
maximum or by nearby gas flowing into or out of that location. If velocity perturbations result 
from spiral structure as described in density wave theory (17), then Fig. 1 will not represent the 
true perturbed surface density but will still be a reasonable representation of the spiral pattern 
©. A dynamical model that self-consistently incorporates density and velocity perturbations is 
required to compensate for this effect, but is beyond the scope of this paper. 

We also map the thickness of the HI gas as a function of position by calculating the second 
moment, T, of the density distribution p for each point in {R, (p) . When we perform the modified 
unsharp masking on T{R, cp) to create r(i?, 0), spiral structure is evident (Fig. 2). The quantity 
r is a direct equivalent to 11 as defined in Eqn. 1 . The value of V ranges in magnitude from 0.32 
to 3.29, although the dynamic range of the majority of the points lie within 0.55 to 1.5. We plot 
the n = 1.1 overdensity contour from the surface density perturbation map (Fig. 1) on top of 
the thickness perturbation map (Fig. 2), showing that there is a good match between the arm 
positions as calculated from the surface density and the thickness; the thickness of the HI layer 
is smaller in the arms than in the rest of the disk. The correlation remains good even if we do 
not apply a thickness filter to the data, or if we use another definition of the thickness 0I77T). 
The two maps are not independent, however; changing the density distribution in the arms will 
necessarily produce changes in the measured thickness unless the perturbations are distributed 
in the same way as the initial distribution. The alignment of overdensities with regions of 
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reduced thickness was suggested previously 01 • This alignment has not been observed in other 
galaxies because surface density maps are most easily made for face-on galaxies where there is 
no information about the thickness of the gas layer. 

The radial profile of the HI disk has been a matter of controversy for many years. A sharp 
falloff in HI emission as a function of velocity has long been known dT^ . but this need not 
correspond to an abrupt radial cutoff in the disk density (179b . Velocity dispersion will cause 
features to be smeared along the line-of-sight by confusing the velocity-distance transformation, 
resulting the radially elongated features near the edges of maps (Fig. 3). The radial extent of 
the spiral arms provides a minimum cutoff radius for the Galactic gas disk; in other words, it is 
not possible for the gas to have spiral structure beyond where the HI disk ends. This radius is 
only a lower limit because it is possible that there is gas beyond where the spiral structure ends 
that does not participate in the spiral structure, or that past some radius the arms are too weak 
to be detected by the unsharp masking. Both the surface density and thickness perturbation 
maps (Fig. 3) change from spiral patterns to features elongated along the line-of-sight near 25 
kpc Galactocentric radius. This is most clearly seen in the south; the transition radius is not 
immediately obvious in the north. Thus, the HI gas disk must extend to at least 25 kpc from 
the Galactic center in the south, about three times the Sun-Galactic center distance. A related 
conclusion is that gas within the cutoff radius is kinematically settled into a disk; otherwise it 
would be unlikely to respond to the spiral density waves. 

It is useful to fit four-armed models to our density perturbation map. We use logarithmic 
spiral arms that start at the Solar circle: 

log(i?/i?o) = - 0o) tanV- (2) 

where ^ is the pitch angle and 0o is the Galactocentric azimuth at the Solar circle. Our fitting 
method is designed to trace the regions of gas overdensity. For each of the four arms apparent 
in Fig. 1 , we investigate an evenly spaced grid of these two free parameters for ranges of values 
that connect the overdense contours. For each combination of ip and (f)o, we linearly interpolate 
the value of 11 for the locus of points along each arm. Any points that fall in the excluded 
regions are ignored. We use the median of the list of interpolated values as a measure of the 
goodness of fit for each curve. In this scheme, arms with values of and cpo that trace overdense 
regions will naturally have a large median, and thus a large goodness-of-fit. The best-fit values 
of tp and 00 for each of the four arms are given (Table 1). Other fits that connect a different 
set of features in the map could be drawn because assigning a unique arm pattern to a map is 
not possible. We find pitch angles for the outer arms in the range 20°-25°; this is larger than 
the value oftp^ 13° averaged over a variety of tracers S2U^ . This does not necessarily imply a 
disagreement, however, as the arms could be unwinding in their outer regions. 

Various models of the locations of the arms have been proposed. We compare our map 
to a model derived from regions of ionized hydrogren {21 22. 23); the model consists of two 
pairs of mirror symmetric arms following logarithmic spirals. We denote this as the symmetric 
model (Fig. 4). The symmetric model fits 11 reasonably well over much of the southern sky; the 
agreement is poor in the north where the spiral structure is less prominent, possibly becayse of 
the larger thickness of the northern gas {11). Gas that is dynamically warmer is less likely to 
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Figure 2: A contour plot of the perturbations in the gas thickness, T{R, (j)). Colored regions 
have reduced thickness compared to the local median, whereas grayscale regions have larger 
thicknesses. The perturbation levels for the different contour levels are given by the colorbar. 
The solid contour marks the line 11 = 1.1, the same as in Fig. 1, to show the alignment of the 
overdense surface densities with the thinner gas regions. 
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Figure 3: A: A contour map of the modified unsharp masked surface density, Il{R, cj)), out to 
40 kpc. B: A contour map of the modified unsharp masked disk thickness, r(_R, (p). The dashed 
circles are lines of constant Galactocentric radius, marked in kpc. 
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Arm ijj (deg) 0o (deg) 



1 


24 


56 
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Table 1: Fits to pitch angles ^ and angles at which the arais cross the Solar circle (f)Q. The arm 
numbers correspond to Fig. 4A. 



respond to spiral density waves, and the azimuthally averaged thickness of the northern gas is 
nearly twice that of the southern gas at R = 20 kpc. 

There are several places where the symmetric model deviates from the data. For example, 
the arm in the north {R ^ 13 kpc) falls in between two of the model's arms; forcing the arms 
to be mirror imaged pairs is too strong a restriction. Features near the excluded regions could 
result from large scale ordered velocity structure that has not been included in our rotation 
model. Elliptical streamlines with m = 2 could cause such an effect {11). Images of other 
galaxies suggest that the spiral arms may bifurcate into spurs in the outer disk. The structure of 
the Perseus and Carina arms past R^ 20 kpc is suggestive of this behavior. 
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Figure 4: A: The same contour plot as in Fig. 1, with a four armed logarithmic spiral fit 
overplotted. The fitting method is described in the text. Other fits that connect different features 
are possible. B: The same contour plot as in Fig. 1, with the four armed symmetric spiral model 
overplotted (23 y The solid lines represent the model over its claimed range of validity; the 
dashed lines are an extension beyond that range. The unlabeled short line near the Sun is the 
local Orion arm. The width of the model arms is arbitrary. 
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